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Craps et al. recently proposed to study a simple cosmological background with a null 
singularity [fl], this cosmology admits a matrix model description, thus lends itself to a 
rigorous study. An earlier example of cosmology with a null singularity is proposed in BB] 
and is later studied by many authors [Af], they find that this singularity is highly unstable 
due to gravitational back-reaction. However, the model of Craps et al. seems to avoid this 
problem. 

The construction of the model in [fl] is simple. One starts in type ITA string theory 
with a flat string metric and a linear dilaton background ¢ = —Qa*, where z is a light- 
like coordinate. The linear dilaton does not require modifying the critical dimension, since 
its linearity is along a null direction. Although the string metric is flat, the Einstein metric 
is nontrivial: 

ds? = e@*"/2(_2dxtda~ + (dz’)?). (1) 


For a positive Q, the metric contracts to a singularity at zt = —oo, this is actually 
a curvature singularity. The corresponding 11 dimensional M theory metric also has a 
singularity at «+ = oo. The authors of [[] show that the singularity at xvt = —oo lies in a 
finite geodesic distance, while the other singularity is at infinite distance. This background 
preserves half of total 32 supersymmetries, one expects that there is a control over the null 
singularity. 

Since the string metric is flat, this time-dependent background admits a matrix string 
description, in which the Yang-Mills coupling constant is time-dependent. In fact, the 
authors of show that this Yang-Mills theory can be regarded as one with a constant 
coupling on a world sheet with a time-dependent metric. 

It goes without saying that this is an important observation, and it may provide the 
first example in which we can study a time-dependent background in a controlled fashion. 
Thus, it is interesting to ask whether this model is unique, or it has many cousins. In this 
note we will see that indeed there is a large class of such models. 


We shall limit ourselves in M theory in this note. Consider metric 
ds? = e?r" (_odrtdx~ + (dx*)?) + ebt" (dg)? (2) 


where i = 2,...10 — d, a = 11 —d,...10, namely, there are d coordinates x°, the total 
dimensions of spacetime is 11. The metric of [I] is given by taking d = 1, a = Q/3 and 
B = —2Q/3. X'° is taken as the M theory circle, thus the metric reduces to the IIA flat 


string metric with a linear dilaton background. 
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To find a solution, we use an orthonormal basis 


“il + wit $ + + 
et =e drt, eb =e% dr!, el =ef? dr’. (3) 


The non-vanishing components of spin connection are 


-azt i 
Wi. = oe e", wip = ae? e, Wap = Be e°. (4) 


The non-vanishing curvature 2-forms are 
2 —2art i 
Rip = ae 7% e et, 


Ra+ = b(2a — Be 202" et Aer. 


The only non-vanishing component of the Ricci tensor is R~+ given by 

R44 = [(9 —d)a? + dB(2a — B)?" (6) 
and the Einstein equation R++ = 0 has two solutions 
Ja. (7) 


For d = 1, choosing the minus sign in (M) reproduces the background considered in [i]. 
We now show that the background given by eqs. (P) and (M) not only is a solution, but 
also preserves half of supersymmetries. The only interesting SUSY transformation is that 
for the gravitino 6V,, = D,e. The components of spin connection of interest are w+, wi 
and wa: . 
w+ =—2a-—2ay-yt, 44 = 2a y", 8) 
8 
Wa = 2Be b- yant, 
Since 
1 =g at a 
[D+ Da] = 5B(B — 2ajeP“" yt, (9) 


The compatibility condition for the constraints D,¢ = 0 is yte = 0. This condition 
eliminates half of the components in e. The constraints D;e = Dac = 0 are solved if € is 
independent of x’ and x°. Finally, D}e€ = 0 is solved if € = exp(5ax™)n for a constant 7. 
We conclude that the unbroken supersymmetry is parametrized by 7 with the constraint 
y*n = 0. Thus, just like the original background of [fl], our more general metric also 
preserves 16 supersymmtries, this should be enough to guarantee that our proposed matrix 


model to be described shortly is a valid description of dynamics over this background. 
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In the case d = 9, 8 = 0, 2a, although [D,, Da] = 0, but [D4, Di] = —faPayiyt, again 
we obtain the same unbroken supersymmetries. 

In fact, (Q) is only a special case of a more general class of solutions preserving half 
of supersymmetries (to be discussed in the end of this note). Before we study the matrix 
model for (P), let us discuss the geometric properties of (P). Apparently, the sign of a and 
its absolute value can be changed by changing xt and x~, so we always assume a > 0. As 
xt — —oo, the factor ezant approaches zero, the transverse dimensions z’ shrink to zero 
size, this is the big bang point with regard to these coordinates. This singularity locates 
a finite distance away in view of the affine parameter, indeed, let dX* = e202" dat the 


first term in the metric (P) becomes 
—2dX*da-, (10) 


thus XT is the affine parameter for the null geodesic x7 =const. The big bang singularity 
occurs at XT = 0, 

The geometry of other transverse dimensions x° depends on the choice (3 in (F). We 
name the choice 3 = (1+3/Vd)a case 1, and the choice 3 = (1—3/V/d)a case 2. For case 
1, 8 > 0, so all z shrink at xt = —oo. For case 2, 3 < 0 except for the d = 9 case, thus 
all x° get to infinity at x = —oo and shrink to zero size at r* = oo. d = 9 is special, in 
this case, all 9 dimensions x* do not evolve in time, and one can redefine «* such that the 
metric is Minkowski. Supersymmetry is also enhanced, since there is no longer constraint 
yte = 0 following from (Q). 

The model studied in [fl] is a special example of case 2 when d = 1. As in [i], we can 
compactify xt? on a circle to obtain type IIA string theory. The string metric and the 


dilaton ¢ are related to the M theory metric through 
ds? = e~ 2/3 qs2, + et? (dx!?)?, (11) 
we obtain 
o= Bet, ds?, = e(20+8)a* [—2datda~ + (dx*)?] + e882" (dy %)?, (12) 
The 10D Einstein metric reads 


ds = e(2at8/4)a* [—2da* da~ + (dz*)?] + e(9B/A)@™ (ga)? (13) 
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Since the string coupling “constant” gs = explo) = exp(2G2*), for case 1, strings are 
weakly coupled at big bang, while all transverse dimensions are zero-sized, so it is not 
clear whether we can trust the perturbative string theory. As we shall see shortly in the 
matrix model, spacetime breaks down and we need to employ a full non-abelian description. 
At later times, strings become strongly coupled, while all dimensions expand, we shall see 
that a rather simple theory emerges, that is, only abelian degrees of freedom survive. 

For case 2, strings are strongly coupled at big bang. transverse dimensions x’ start 
with a zero size, while transverse dimensions z“ have infinite size and contract as time 
evolves. This picture in valid both in terms of the string metric as well as the 10D Einstein 
metric. We will study in more detail the spacetime properties in the matrix model later. 

The string frame metric ([2) is in general non-flat. Ø never vanishes, so to get a 
flat metric, d = 1 is necessary, and in this case one chooses @ = —2a, this is the special 
case considered in [fl]. More generally, one obtains a world-sheet theory with an action 
explicitly depending on time. One can attempt to quantize the string in the light-cone 
gauge, to do so, introduce the new light-cone coordinate dy* = exp((2a + 3)a*)dx*, the 
first term in the string frame metric ([2) becomes —2dytdx~. The lignt-cone momentum 
p_ = —p* conjugate to x7 is conserved. In the light-cone quantization, we can set the 


light-cone gauge yt = 7, the bosonic part of the world-sheet action reads 


1 

g= f drda[e@r+)2" 9 X Ia X? + e998" A°X%9,X 4, (14) 
Ara! 

where the period of o is 2ra'p™, and exp((2a+ 8)xt) = (2a+)r. As xt — —oœ, T > 0. 

Since 7 starts at a finite point, it is more useful to use the old light-cone coordinate xt = t, 


the action becomes 


1 
SS f Gide O" X40, ke AX]: (15) 
TQ 


The time-dependent coefficients may be interpreted as the effective tensions. For transverse 
coordinates XŻ, the effective tension is T; = sire 4028 )t for coordinates X“, the effective 
tension is Ty = Poi: mle In case 1, both effective tensions get to zero at big bang, 
while the string coupling also gets to zero. Since the string spectrum becomes very dense, 
the usual free string picture does not apply. In later times, effective string tensions get 
large, and the string coupling constant also becomes strong, we will have to deal with a 
strongly coupled massless sector. In case 2, 4a + 23 > 0 except for d = 1, our above 


analysis still applies to T;. 2a + 48 = 6a(1 — 2/Vd), for d > 4, Ta becomes small in 
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earlier time, the same as case 1. For d < 4, Ta becomes large in earlier time, thus these 
transeverse dimensions are effectively frozen. d = 4 is a special case, T, is independent of 
time. 

The vertex operator of a string state in the background (IJ) in general is quite com- 


plicated. For instance, consider a massless scalar satisfying equation of motion 
On (e~??./—ggt”0,®) =, (16) 


As function of «+, the scaling of components g% is different from the scaling of components 
g*~ and gř, so there is no simple plane wave solution in general, which implies that the 
vertex operator of this massless scalar field is not simple. However, we can consider special 
cases when ® is independent of x“, in this case, a vertex operator V = exp(ikyat + 


ik_x~ + ik,x*) must satisfy the on-shell condition: 


k-(2k} —iy)-k? =0, y=(9-da+dZ, (17) 
which require a imaginary part of k}: Sk_ = 7/2, thus in the vertex operator, there is an 
exponential factor 

enn", (18) 


Since y = (9+ 3Vd)a > 0, this exponential factor always blows up at xt = —oo for d < 9. 
The effective string coupling constant for these states is geff = gs exp(—3727). Dimen- 
sions x° for these states are effectively compactified, thus when we discuss interactions 
among these states, the space-time dimensionality is 11 — d, and there is an effective string 
tension in action (1), or al pp = a’ exp(—(4a + 26)x*). The effective Newton constant is 


9—d)/2 


Geff = G2 pou’ , as a function of x*, it scales 


Geff ~ ete —(9-d)(2a+f)2* (19) 


Now, y is positive and (9 — d)(2a + 8) is non-negative if d < 9, the effective Newton 


+ = —oo, the perturbative sstring picture is not valid at least for 


constant blows up at x 
these states. 
Next, we study the matrix model. It is not necessary to compactify any dimension of 


(P). In the usual flat background, the Matrix Theory action reads [f] 


1 : : ; . 
S= far (KO + FIX, XI)? + 107 Di0 — RO? yX, a) , (20) 
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where R is the longitudinal cut-off, or it may be viewed as the radius of x~ in the DLCQ 
M theory [$]. For simplicity, we set the M theory Planck length lp = 1. A derivation of this 
matrix action is given in [P] (see also [[0] ). With our metric (Ø), it is straightforward to 
write the corresponding matrix action. To do this, we need to use the light-cone coordinate 
in (10), and identify X* with 7, the world-line time of D0-branes. In the matrix model, 
we apply the M theory metric (Q) directly. Since the action is rather lengthy, we separate 


the action into the bosonic part and the fermionic part. The bosonic part reads 


1 1 R 
g / dTa ; eart (D, xY ; e2bet (D, xa)? ic etar X, X1] 
(21) 


R R , 
A mesa w eba X", X lla $ oo Foe aig? 


Note that xt appearing in the above action is the old light-cone coordinate. The fermionic 


part reads 
Sp = f dr{i0™ D-0 — Re? 67 yi[ Xt, 0] — Re?” 0Tya| X", 6] }. (22) 


It is rather awkward to use 7 as time, since it has a finite beginning. Let us switch 


back to the coordinate x™ and on the world-line identify t = x*, thus dr = exp(2at)dt, 


we have 
1 i\2 1 2(B—a)t ava E 6atryi 174 
Sp = | dtTr{—(D_X")* + —e (DiX) + =e X", X?] 
2R 2R 4 (23) 
+ E a+4a)t ya xe B (4a+29)tj yi xh 
4 X 9 7 3 
and 
Sr = I dt{i0T D,@ — Re®“*97 y,[X*, 0] — ReCot tot y |X", 0]}. (24) 


Before study the simplest properties of this matrix model, let us show that we can 
recover the matrix string action of [I]. In this case, there is only one X%, call it X'°. 
The compactification scheme is given in |I|. Up a dimensionful parameter, we replace 
the trace in the matrix action by f doTr, the commutator R[X!°, X"] is replaced by the 


covariant derivative iD, Xt. Finally, use 8 = —2a, we find 
Sp= dtdoTr{ D Ky ke 2 p2 pyar XIP) (25) 
B 2R a IRIs to 478 ’ ’ 
and 
Sr = I dtdaTr[07 o" Da0 — Rgz 10" 4;[X', 0], (26) 
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where gs is the time-dependent string coupling gs = exp(—3at). Our action is identical 
to that in up to an identification of a dimensionful parameter ls. This matrix string 
theory can be regarded as a 2D Yang-Mills theory with a time-dependent coupling or a 
2D Yang-Mills theory with a constant coupling in the world-sheet metric 


ds? = g,°(—dt? + do’). (27) 


To study the properties of matrix model defined by (23) and (24), we consider cases 
1 and 2 separately. 

e Case 1, 8 = (1+ J)a. 

The kinetic term of X’ is always simple, as is the kinetic term of 0. The coefficients of 
the remaining terms all vanish in the limit t — —oo, so there is no constraint arising from 
these terms, this implies that all matrices are fully non-abelian. On the other hand, as t — 
oo, these coefficients blow up, thus, all matrices must commute with one another, the only 
surviving degrees of freedom are diagonal elements. Moreover, X* must be independent of 
time in this limit, so X° become frozen abelian moduli. Recall that in the string picture, 
we found previously that the effective tension T, becomes infinitely heavy, even heavier 
than T; for large t, this is related to the fact that X“ become moduli in the matrix model, 
while XŻ still have dynamics. If we compactify some of the transverse dimensions, the 
story becomes slightly more involved. For instance, compactifying X!° on a circle, all 
other matrices become function on a circle ø. In the limit t — oo, although they have to 
be diagonal, they are not always periodic functions of a, the eigen-values can get permuted 
after circling along g, these twisted sectors describe strings of various lengths. 

e Case 2, 8 = (1 — qe 

6a and 4a+2 are always positive, at big bang, there is no constraint on commutators 
[X*, X1], [X*, X%], [X*, 6] and [X“%, 0]. On the other hand, as t > œo, these commutators 
are forced to vanish. 2(3— a) is always negative, so at bing bang, X° must be independent 
of time, they become non-abelian moduli in the model. 2a — 43 = 6a(1—2/Vd), d = 4 
becomes the critical dimension. For a larger d, there is no constraint on the commutators 
[X*, X°] at big bang, and they have to vanish as t — oo, thus, for d > 4, all matrices 
commute in this limit and only abelian degrees of freedom survive. For d < 4, [X%, X°| 
have to vanish at big bang, together with the fact that X° are independent of time, these 
matrices become abelian moduli at big bang, this fact is also reflected in our previous 


analysis in the string picture, where we found Ta become infnitely heavy at big bang. 


T 


When d = 4, X*% remain nonabelian all the time. In the string picture, recall that when 
d = 4, the effective tension Ta is constant. 

It is also of interest to study compactification. Just as in the flat background, we do not 
know how to write down matrix model action if we compactify more than 5 dimensions. 
We may choose to compactify all X°, or some of X%, or some of X% and some of X’. 
For simplicity, let us compactify all X° on a torus T¢. Replacing R[X*%, X*] by iD, X’, 
R?[X°, X4| by —Fup etc, the bosonic part of the matrix model action reads 


1 : 1 : 1 
— d+1 D xX? 2 (4a+28)t D x? 2 2(8—2a)t p2 
oe [« igp A =p (DA taga" a 
1 R ae 
ae i ee 4 zo" [X*, X1. 


The analysis of this action is similar to our previous analysis of the action without com- 
pactification, a statement there can be simply translated to one for the action (R8). For 
instance, demanding | X°, X?] = 0 is translated to Fap = 0, that is, the spatial connection 
must be flat. 

It is interesting that the previous noticed “critical dimension” d = 4 corresponds to the 
situation that the action (28) is no longer complete. Compactification on a torus T4, the 
complete matrix model is to be given by the world volume theory of coincident M5-branes 
[2]. Of course, that theory will also be time-dependent, and hopefully one can figure out 
more details of the time dependence by studying field theory behavior of (28). Likewise, 
action (R8) is again incomplete for compactification on T 5 in this case the complete theory 
is given by the little string theory [13]. 

Unlike the case when d = 1, the Yang-Mills theory (28) can not be simply interpreted 
as a theory with a constant coupling on a world volume of a nontrivial metric. To interpret 
action (28) as Yang-Mills theory, we need to introduce both a nontrivial metric on the world 
volume as well as a time-dependent Yang-Mills coupling. Since a factor Iya appears in 
the coefficient of F2, and a factor gẹ appears in the coefficient of [X+, X1]?, we find for 


d>1 


ds%yy — eT (4a+26)t 42 te EFT 40428)t (gga)? 
29 
a = e7 T% (4a+28)t+2(a—8)t ey 


where the subscript WV stands for world volume. the world volume metric components 
Jaa Must be all same, there are only 3 independent functions in the geometry and the 
Yang-Mills coupling, however, we need to match 5 different kinds of coefficients in (28), 


so it is nontrivial that there is a solution as presented in (29). This is quite important 
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especially for d > 3, since we need to complete the action (28) by introducing more degrees 
of freedom, this is possible only when (28) admits Yang-Mills theory interpretation. One 
can check that the identification (29) also works for the fermionic part of the matrix action. 

Finally, we discuss generalizations of background (P). There are two directions to 


generalize (P). The first is to consider 


ds? = —2620?" datdx~ + > e?Pizt (dat)? (30) 


The Einstein equations are solved provided 
a? + X Bila — Bi) =0. (31) 


This background again preserves 16 supersymmetries of the form e€ = exp(saxt)n, with a 
constant 7 satisfying y*n = 0. 


The second direction to generalize (Q) is to consider metric of the form 
ds? = ef") (_odat da~ + (da*)*) + eI) (da®)?, (32) 


The non-vanishing curvature 2 forms are 


R; = (F°? _ fe? e A et, 


(33) 
Ra = (2f'g' _ g? P ghe? es IN et. 
The Einstein equations are solved if 
(9 — d)(f? — f") + d(2f’g’ — 9? — 9") =0. (34) 


Once again, there are 16 unbroken supersymmetries, parameterized by € = exp(5 fn, 
yin =0. 

Needless to say that equation (B4) admits infinite many solutions. Two special solu- 
tions deserve attention. One solution is when g = 0, in this case f’? — f” = 0. All the 
curvature 2 forms vanish thus it may appear that we obtain a flat space solution. This is 
almost case except for a simple twist. The solution of f’? — f” = 0 is exp(2f) = 1/(x*)?, 


by changing coordinates dyt = dx*/(x*)?, y~ = x7 we obtain a metric 


ds? = —2dyt dy” + (y*)?(dz")? + (dx). (35) 
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If one of x’ is periodic, we find a nontrivial background similar to the orbifold discussed 
in ØB]. Of course if all zê are noncompact, (BJ) is equivalent to the Minkowski space. 
Similarly, one may set f = 0 and find exp(2g) = (xt)?, a nontrivial orbifold results 
provided one of xz“ is compactified. One may attempt to construct a matrix model for 
the metric (B5) too. However, time r = y* seems to have to terminate at r = 0. One 
may try to avoid this problem by going back to the original light-cone coordinate x* and 
take t = xt in the matrix model, although this helps to eliminate the coefficient T° in the 
kinetic term of X* but it also introduces a coefficient t? in the kinetic term of X°, again t 
must be terminated at t = 0. Also the supersymmetry parameter € = Vx+n does not exist 
beyond xt = 0 since € has to be real. This problem may cause disease in the definition of 
the matrix model, and may be a reflection of the gravitational instability caused by a test 
particle. 
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